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We study the evolution of the quantum state of n photons in m different modes when they go 
through a lossless linear optical system. We show that there are quantum evolution operators U 
that cannot be built with linear optics alone unless the number of photons or the number of modes 
is equal to one. The evolution for single photons can be controlled with the known realization of 
any unitary proved by Reck, Zeilinger, Bernstein and Bertani. The evolution for a single mode 
corresponds to the trivial evolution in a phase shifter. We analyze these two cases and prove that 
any other combination of the number of photons and modes produces a Hilbert state too large for 
the linear optics system to give any desired evolution. 
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I. QUANTUM OPTICS IN 
PHOTON-PRESERVING LINEAR SYSTEMS 

There are many optical elements that can affect the 
quantum state of light. Elements that preserve the num¬ 
ber of photons are particularly interesting in quantum 

f tics and in applications to quantum information [ll- 
Linear, lossless, passive systems have received a great 
deal of attention since the demonstration that, combined 
with measurement, they can be used to build a univer¬ 
sal quantum computer [^. Recently there has been a 
revived interest kindled by the result that the output 
statistics of linear optics multiports cannot be accurately 
predicted in a classical computer efficiently unless sev¬ 
eral well-founded computational complexity hypothesis 
are false 

In this paper, we study the behaviour of optical sys¬ 
tems that act on n photons in m different modes. We 
call m X m multiports to the optical systems of interest. 
The evolution of the state of the photons can be charac¬ 
terized from the scattering matrices S used to describe 
TO-ports in classical electromagnetism. We stick to the 
port denomination for the intuitive picture it gives, but 
the photons can really be in different orthogonal modes. 
The key is that two photons in different modes are per¬ 
fectly distinguishable and do not interfere. The simplest 
example is a system with photons travelling in different 
paths, but we can also imagine photons in orthogonal po¬ 
larization states or which have orthogonal orbital angular 
momentum states. 

The inputs to our system are a combination of states 
with rii photons in a mode with index i, denoted by |ni)j. 
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Eor a system with a total number of photons n, all the 
possible input states can be described as a linear combi¬ 
nation of states 

l'^) = \ni)i\n2)2---\nm)ra ( 1 ) 

with ni -|- n2 -l- ... -|- Um = n. Linear optics multiports 
present at their output a linear combination of states of 
the same form. 

The evolution of a photonic quantum state in our sys¬ 
tem can be specified from a unitary matrix U so that 
l^out) = R|dfin)- The classical scattering matrix S is 
enough to characterize the evolution of any number of 
photons entering the multiport. Both S and U must be 
unitary matrices as they describe systems that conserve 
energy and the total probability, respectively. 

The step from S' to 17 depends on the number of 
photons. If we take the basis composed of the num¬ 
ber states of Eq. o, the element of U that describes 
the transition from I'l'in) = |u2)2 ■ • ■ 

I'kout) = 1742)2 ■ ■ ■ Wm)m determined from 

Kli («2l2 ■ • ■ (^mL U 1^1)1 I’^2)2 ■ ■ ■ \nm)m^ which has a 
value 

Per(»S^ijq out) 

• n^! • • • n'^\ • ni! • 742! • • -n™! 

In Eq. (m, Per(Sin,out) is the permanent of a matrix 
Sin,out with elements Si^j from S such that each row in¬ 
dex i appears exactly 74' times and each column index j 
is repeated exactly rij times [ 1 , . 

Alternatively, we can write our number states from 
their creation operators so that \ni)^ = | 0 )j and see 

how the operators transform. For a linear optics multi- 
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port, we know Q the creation operator a| evolves into 

m 

( 3 ) 

i=i 


The size of the scattering matrix is a function of the 
number of inputs and outputs of the optical system. S is 
an m X TO matrix, whereas U is an M x M matrix, with M 
the size of the Hilbert space that contains all the possible 
configurations of n photons divided into to modes. These 
different states form a complete basis of the state space 
and their number is equivalent to the number of ways of 
placing n indistinct balls in to different boxes, which is 
the combinatorial number 


M = 


/to + n — 1 

V n 


(to + n — 1)! 

(to — 1)! n! 


(4) 


We can generate all the possible states recursively if we 
assign a photon number i from 0 to n to the first mode 
and then generate all the possible states for the n — i 
remaining photons in the rest of the modes. By the time 
we arrive to the last mode the assignment is trivial and 
we can repeat the procedure until we have a complete 
list. 


II. UNIVERSAL QUANTUM 
TRANSFORMATIONS 

We say we have universality if, for our number of pho¬ 
tons n and the number of modes to of our system, we can 
generate any desired quantum evolution U in the state 
space of all the possible distributions of the n photons in 
the TO modes. 

In this paper, we show there are limitations to the 
quantum transformations U we can create from a linear 
optics multiport. While we can implement any desired 
unitary scattering matrix S using only beam splitters and 
phase shifters Q, a tailored S can only produce any ar¬ 
bitrary U in a limited set of cases. 

This is a problem different from finding a universal set 
of gates for quantum computation. In most linear optics 
implementations of quantum computing we restrict our¬ 
selves to only a subset of all the possible quantum states 
and there is some kind of postselection. 


III. DEGREES OF FREEDOM AND 
UNIVERSALITY 

The main result of the paper is a proof that there is a 
necessary condition for universality which is only satisfied 
in a limited number of cases for which there are explicit 
ways to describe how we can generate any desired U. 

The basic argument is that the degrees of freedom we 
have when we build the multiport must be at least equal 
to the degrees of freedom in the photonic Hilbert space. 


Otherwise, there will be transformations that are impos¬ 
sible to perform. 

Lemma 1 A linear optics multiport with m inputs can¬ 
not be used to give all the possible quantum evolutions in 
the state space of n photons in to distinct modes unless 
TO > M, where M is the dimension of the Hilbert space 
of the photonic states. 

Proof The unitary group Uirn^) contains the mxm 
matrices S that describe the linear optics system and 
the unitary group U{M'^) contains the M x M matrices 
U that describe the quantum evolution of the photons’ 
state. Using the expression of Eq. we can define an 
homomorphism ip : S ^ U which maps U(rn?) to U{M‘^) 
We can only reach all the matrices in if ip 

is surjective, which for our unitary groups is equivalent 
to ask for p to be an epimorphism. The homomorphism 
can only be surjective if the dimension of the domain of 
p is at least as large as its codomain. In our problem, 
the condition is m? > M^, which, for the ranges we are 
interested in, reduces to the necessary condition for uni¬ 
versality 


TO > 


TO -|- n — 1 
n 


(5) 


The intuition behind this result is that we have only a 
limited number of degrees of freedom when we build the 
linear optics system. If the target state space is too big, 
we cannot reach all the possible matrices U. 

In the following sections, we show that in all the cases 
where necessary condition is met (n = 0, n = 1 and 
TO = 1), there is also an explicit way to find any desired 
unitary. For n > 1 and to > 1 we prove it is impossible 
to implement all possible unitary matrices U using linear 
optics alone. 


A. The vacuum state is always taken to the vacuum 


The first trivial result is that linear optics preserves 
the vacuum state with zero photons. This is obvious as 
a passive linear optics multiport cannot create photons, 
but can also be deduced from the necessary condition of 
Eq. ([5]). Our Hilbert space has a dimension 


/m-\- n — 1\ 
\ n J 



( 6 ) 


and TO > 1 for any linear optics system, which will have, 
at least, one input. There can be many unused degrees of 
freedom. With no photons the exact configuration of the 
linear optics multiport is irrelevant and we can choose 
different scattering matrices. 
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B. Systems with one port are equivalent to a phase 
shifter and trivially give universality for any number 
of photons 


When m = 1 we have a similar situation. For any 
number of photons n 


M = 


m + n — 1 
n 



(7) 


and the necessary condition of Eq. (ED is met with m = 
1 = M. The interpretation is also clear. If we have 
only one mode, the only allowed physical operation is a 
phase shift which is equivalent to a 1 x 1 unitary matrix 
S = whose only element is a root of unity. The 
linear optics system can only be a phase shifter. The 
evolution for n photons is then a phase term e™'^. We 
can use our degree of freedom (j) to give any output phase 
shift we want and we have universality. 

However, in a quantum state we cannot observe a 
global phase shift. Phase can only be determined when 
compared to a reference, like in interference between 
states. This is similar to the definition of voltage, where 
only differences of voltage have a physical meaning. The 
output state |n)^ is equivalent to the input state 
\n)-^. There is no measurement that can distinguish 
between these two states. For a set of measurement 
operators {Eq}, we obtain outcome o with probability 
p{o) = {n\^ElEo I n)-^ for |u)]^, which is the same result 
we get for (u|^ e“®‘^if]ifoe®‘^ | n)^. In this case, any uni¬ 
tary matrix U is, really, equivalent to the identity matrix 
L 

Notice that the equivalence disappears if we have a 
reference state. If we had one photon in a reference path, 
the effect of the phase shifter could be observed with a 
well designed measurement. But then we would be in a 
different case with m > 1. 


C. Linear optics multiports can give any possible 
quantum transformation for one input photon 


beam splitters and phase shifters or only one kind of 
beam splitter for m > 3 I,®- 

This result cannot be used to build a scalable univer¬ 
sal quantum computer. If we want to implement an al¬ 
gorithm acting on q qubits, we need to use 2'J paths to 
generate all the possible states. This exponential growth 
prevents a generalized use of linear optics with one pho¬ 
ton for quantum computation. 


D. Linear optics alone cannot give any desired 
quantum transformation for more than one input 
photon in more than one mode 


Apart from the limited results of the previous sections, 
in general, linear optics multiports cannot give any de¬ 
sired unitary evolution. 

Theorem 1 A linear optics multiport with m > 1 inputs 
cannot he used to give all the possible quantum evolutions 
in the state space of n > 1 photons in m distinct modes. 


Proof We consider all the cases with m > 1 and n > 1. 
From Eq. o 


M = 


TO -I- n — 1\ (to -|- n — 1) • • • (to -I-1) • TO 
n J nl 

= ”- rim -■ ® 


We can write the dimension of the photons’ Hilbert space 
as 


M = m 


TO-l-1 TO TO-|-n—1 


n 


m — 1 k 


( 10 ) 

with a product of terms 1 -I- > 1 if to > 1. For u > 1 

and to > 1 there is at least one such term in the product 
and it is immediate to prove M > m which violates the 
necessary condition of Eq. ED- ■ 


The next interesting case is the evolution of a single 
photon in an arbitrary multiport with to > 1 ports. Here 

frn + n — 1\ (m\ 

M=( ^ j = (i)="' («) 

and we fulfill the necessary condition with to = M. One 
basis of the Hilbert space of the photon states is the basis 
of elements |f) where i = 1,...,to is the index of the 
mode our only photon is in. State |f) corresponds to a 
column vector filled with zeros and a 1 entry in row i. 
From the definition in Eq. (ED we can see that, for one 
photon, (j\U\i) = Sj^i- The permanent Per(5'in_out) is 
exactly the only element of the matrix iSin.out = 
and U = S. 

We can then implement any desired unitary directly 
by choosing the appropriate matrix S. There are con¬ 
structive methods to implement any unitary S using only 


IV. COMMENTS AND EXAMPLES 

We have shown that, except for a few restricted cases, 
linear optical systems cannot be used to give any desired 
quantum evolution for n photons divided into to optical 
modes. We have given a necessary condition for uni¬ 
versality and proved that when the condition is satisfied 
there are explicit constructions for any unitary evolution 
U we require. 

It is still open how severe this restriction is. In the 
condition m> M, the growth of M as suggests 

a smaller number of achievable operators U for higher 
values of n and to. However, there can be important 
limitations even for small state spaces. We can show 
that in a simple example with two input ports and two 
photons. 
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The linear optics system is determined by the unitary 
matrix 


input. If we use the general expression in Eq. m , we 
can see the output state for an input |11) is 


From Eq. we can see an input state 1 ^ 2)2 has 
an output 

^i==(5nal + 52ia^)'‘n5i2al + ^ 220 ^)-= 10)^ lO)^ . 
\'ni\n2l 

( 12 ) 

For the basis {|2)^ [ 0)2 , |0);^ |2)2 , |l)i 11 ) 2 } and defining 

oj, |0)i|2)2 = |02)= 1^1 j , 

/o 

|1)J1)2 = |11)= I 0 

the unitary matrix that gives the evolution of the pho¬ 
tons’ state is 


/ V 2 S 11 S 12 \ 

l^out) = 72^21^22 • (17) 

\ 5 'i 1 <S '22 + > 5 'i 2 *S '21 / 

Imagine we want to obtain the output state |20). We 
know this is impossible because it would require the ma¬ 
trix of Eq. (HU), up to a global phase. We can, instead, 
search for the closest possible state, as measured from 
the overlap 

|(20|t/|ll)|2 =2|,Snn,Si2|' = 2|5ii|2(l-|5n|'), (18) 

where we use S is unitary and therefore |5'iip -I- |/S'i 2 p = 
1. We would like to get lyes'll5121 = 1, but we must 
settle with maximizing 2|5iip(l — |5iip). The entry is 
maximized for |5iip = 1 with a maximum overlap 1, 
which is exactly the case in the Hong-Ou-Mandel exper¬ 
iment. This example shows the limitations can be severe 
even for values of M slightly above m, like our example 
with M = 3 and n = 2. 

We can also wonder if the results are valid outside Fock 
states. However, other states, like coherent states 




U = 


o2 

>->21 

V25ii52i 


5?2 V25ii5i2 

5I2 V2S21S22 

\f2S\2S22 5ii522 + 5i2521 


(13) 


We can review many interesting known phenomena 
from this description. Take for instance the Hadamard 
matrix that corresponds to a balanced beam splitter 


S = 


1 

71 



(14) 


For an input state |11), if we substitute the relevant terms 
in Eq. (IT^ and operate, we get the evolution 


Hill) 


| 20 ) - | 02 ) 
72 


(15) 


This is the simplest example of quantum interference be¬ 
tween indistinguishable photons and it is described in the 
famous Hong-Ou-Mandel experiment [Tl| . 

We can also use this simple example to show there are 
forbidden operations. For instance, the evolution 


U = 


0 0 1 \ 
0 10 
1 0 oy 


(16) 


r)=, 




00 h. 

a 


k=0 ^ 


(19) 


can always be written as a linear superposition of num¬ 
ber states. Linear optics preserves the number of photons 
and we can study separately the evolution for different 
photon numbers. For most of the terms in the superpo¬ 
sition we cannot achieve any arbitrary evolution. Unless 
we only have superpositions of states for which universal 
evolution is possible there will be forbidden operations. 

The restrictions of the achievable evolutions U does 
not mean we cannot produce any desired output state. 
We can always introduce the desired state Ilk) in a lin¬ 
ear system with a scattering matrix S and measure the 
output 1$). The inverse system, with matrix , will pro¬ 
duce an output |4t) for an input |$). Trivially, if 5 = / 
we can generate any output by choosing that state at the 
input. This only shows that arbitrary state preparation 
is equivalent to preparing a known state and being able 
to perform an arbitrary evolution. 

The interest of the presented result lies in the realiza¬ 
tion that certain states cannot be achieved from certain 
inputs. Determining which states can be reached for any 
given input state is left as an open problem that will 
require different methods that the ones presented here. 
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